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Abstract 

An adiabatic approximation is used to derive the binding potential 
between two heavy-light mesons in quenched SU(2)-colour lattice QCD. 
Analysis of the meson-meson system shows that the potential is attractive 
at short- and medium-range. The numerical data is consistent with the 
Yukawa model of pion exchange. 

1 Introduction 

Quantum Chromodynamics is widely accepted as the correct theory of quark- 
quark interactions. QCD effective theories such as Chiral Perturbation Theory 
have successfully described a wide range of hadronic phenomena. Nucleon- 
nucleon interactions, for example, have traditionally been modeled by effective 
meson-exchange theories. Such theories are long-range approximations to QCD 
interactions mediated by gluons and quark-pairs. A derivation of the strong 
nuclear force from full QCD is obviously desirable. 

Lattice QCD provides a first-principles framework for such a derivation. 
However, the nucleon-nucleon system, containing six light quarks, is well be- 
yond the current limits of computation. In fact, at first sight such a computation 
seems impossible — Lattice QCD is providing physical hadronic data consistent 
with experiment at the ~ 5% level. The nuclear physics problem, the inter- 
action between two colour-neutral objects, is an effect two or three orders of 
magnitude smaller. Nevertheless, nuclear models involving static quark clus- 
ters in unquenched QCD jlj and light-quark hadrons in confining QED2+1 
have been examined, and both give encouraging results. An unquenched lattice 
QCD simulation of the nucleon-nucleon system is doubtless some years away. 
In this paper we examine the same fundamental problem in a simpler setting, 
by examining the interactions between two heavy- light mesons ||. 

An elementary problem in molecular physics is the derivation of the bind- 
ing energy of the H2 molecule. A good first approximation can be found by 
using an adiabatic approximation, in which the (slow) nuclear kinematics are 
removed by treating each nucleus as a static force centre, and the remaining 
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(fast) electronic problem is solved as a function of the nuclear separation. This 
generates an effective potential for the two nuclei, which can be used as input 
for the internucleon wavefunctions and energy levels. 

We wish to consider an analogous problem — a 'hadronic molecule' consisting 
of two heavy-light mesons, in which the heavy quarks are treated as static colour 
sources, playing the role of the (slow) atomic nuclei. The gluons and light-quarks 
play the role of the fast degrees of freedom. In addition to the static heavy- 
quark gluon-exchange interaction, this calculation will include the effects of 
interactions between gluons and light-quarks, as well as light-quark exchange. 
The total energy of the system can be found from the asymptotic long-time 
behaviour of the meson-pair propagator. An effective potential is extracted as 
a function of the meson separation. 

The system described above contains only two light valence quarks. As 
such, we arc still quite removed from a direct simulation of the nuclcon problem. 
However, as all nuclei arc hadronic molecules, our qualitative conclusions should 
be somewhat universal. Simulations of MM- and MM-systems are of further- 
interest because of their immediate application to KK phenomena. Studies of 
multi-quark states indicate that the only likely bound four-quark systems are 
mesonic molecular states Q. Two exotic particles, the do (980) and Jo (975), are 
thought to be lightly-bound KK molecules [B . 

This paper describes a lattice simulation of heavy-light mesonic molecules, 
using SU(2)-colour quenched QCD for computational simplicity. We consider 
both meson-meson and meson-antimeson systems. The next section describes 
the meson operators used in this simulation, and outlines the various contribu- 
tions to the meson-pair propagators. In Section III we list the Lattice actions 
used and the values of simulation parameters, and our results are given in Sec- 
tion IV. 

2 Operators and propagators 

We used a local pseudoscalar meson operator, 

M(x,t) = i)i(x,t)j 5 tp h (x,t) , (1) 

constructed from heavy (h) and light (I) Wilson quarks. With quark propagators 
given by 

G(*,t;0,0) = (0M*,t)^(0,0)|0), (2) 
the zero-momentum meson propagator is then 

G M (t) = Tr t; 0,0) 75*7,(0,0; 2, tbfc 

= T^G h (2,t;0,0)G J t (2,t;0,0), (3) 

X 



2 



The hopping-parameter expansion of the quark propagator gives, in the 
limit of infinite quark mass, 



G h {x,t; 0,0) ocSl >0 JJ t7 4 (0,- 



(4) 



That is, the heavy-quark propagator becomes simply a string of gauge-field 
links in the time direction, and is calculated in one sweep of the lattice. Figure 
1 shows a contribution to the heavy-light meson propagator. 

An operator for a heavy-light meson-pair with spatial separation R is 

M M{R,t) = M{3,t)M{x + R,t). (5) 

The propagator for this system is then 

G M M(t,R) = G D -G E , (6) 

where a sum over all contractions of the quark fields was performed, giving 



G D {t,R) 



G E {t,R) = 



Tr 
xTr 
Tr 



G h (0,t;0,0)G}(0,t;0,0) 
r G h {R,t;R,0) Gl(R,t;R,0) 
G h (0,t;0,0) G\(R,t;0,0) 



G h (R,t]R,0) G](0,t;R,0) 



(7) 



Note that, with the heavy-quark approximation (^J), the separation R is well- 
defined throughout the propagation. Contributions to Go and Ge, the direct 
and exchange terms, are depicted in Figure 2. 

The meson- ant imeson propagator has a similar form, 



Gmm^i R) = Gd - Ge , 



(8) 



where 



Go(t,R) = 



G E (t,R) 



Tr 
xTr 
Tr 



G h (0,f,0,0) G\(0,f,0,0) 
r Gi(R,t;R,0)Gl(R,t;R,0) 
G h (0,t;0,0)G}(R,t;0,0) 



G{(R,t; R, 0) Gz(0,t; R,0) 



(9) 



Contributions to these terms are shown in Figure 3. 

The ground-state energy of the two-meson system is extracted from the 
long-time evolution of the system's propagator |6j. Expanding the propagator 
as a sum over eigenstates of the Hamiltonian |7V) with eigenvalues En, the 
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propagator becomes proportional to the exponential of the ground-state energy 
Eq for long times, 



G(R,t) = (0\d M M(R,t)O M M(R,0)\0) 

= (°\ & mm (£, *) \N) {N | Omm(R, 0) |0) 



JV 

J2\(°\£>mm(R, 0)1 TV) 12 

N 

-E a (R)t 



-E N (R)t 



~ coe^ ^* as t^oo. (10) 

The binding energy V{R) is then defined as the difference between the ground- 
state energy and the masses of the two independent mesons, 

E (R) = 2m M + V(R). (11) 



3 Simulations 

A model of quenched QCD with S'/7(2)-colour was used to simplify analysis. 
The standard Wilson action for the pure gauge field is , 



S G = -/3£(l- iTrUp) . (12) 



Here (3 is the inverse of the square of the coupling strength, and the sum is over 
all plaquettes on the lattice. The Wilson action for light quarks has the form 

Sq^^lpn K nm 1p m , (13) 

with the action matrix given by 

Knm = S n , m ~ K ^ [(1 ~ 7m) UfJ.( n ) S n +fi,m 

+ (l+7^)^(n-A)<5 )■ (14) 

The light-quark propagator is then the inverse of the action matrix, 

G l (n,m)=K^ l (15) 

Simulations were performed on an 8 4 lattice with periodic boundary con- 
ditions. An ensemble of 400 gauge field configurations was generated with a 
Metropolis algorithm at = 2.3. Using the standard prescription for translat- 
ing physical SU(3) QCD to unphysical SU(2) QCD, we used the pion and rho 
meson masses to set the lattice scale M. We found the lattice spacing to be 
roughly 0.2 fm, and the lattice is then ~ 1.6 fm along a side. One expects the 
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radius of the heavy-light meson to be around half that of the pion, or ~ 0.5 fm, 
and so our lattice volume should be large enough to contain the two mesons. 

To ensure statistical independence, 200 update sweeps were performed be- 
tween stored configurations. Quark propagator matrix inversions ( |i"5| ) were per- 
formed with a conjugate gradient algorithm, at a value of k — 0.165 for the 
hopping parameter. 



4 Results 

All propagators under consideration are symmetric in time, since the periodic 
boundary conditions allow contributions from quarks and antiquarks winding 
around the lattice. The mass of the heavy-light meson was taken from fits to 
the meson propagator (g), 

G M (t) = c M cosh (m M (t - T/2)) , (16) 

where T = 8 is the period of the lattice in the time direction. While this is 
an asymptotic relation, valid for long times, in practice fits with acceptable % 2 
were obtained for 2 < t < 6. The average mass obtained from these fits was 

m M = 0.969 ± 0.008 (17) 

in lattice units. 



4.1 The MM-system 

The meson- pair propagator, @, was constructed for values of the separation R 
from to 7 lattice spacings. Given that statistical fluctuations in the individual 
meson propagators will likely be highly correlated to fluctuations in the MM- 
system propagator, we were able to reduce statistical errors by analysing the 
ratio 

C{R,t) = ^M. (18) 

For each value of the separation R, the configuration average of C(R, t) was 
fit with a ratio of cosh functions, 

C( R ,t) = c M ^<^"- V ^ r \ (19) 
cosh (tum t) 

with 

r = t - \ , (20) 

and using the meson mass given in @. The results of the fit for R = and 
R = 1 are shown in Figure 4. 
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The binding potential extracted from ( |ilf ) is shown in Figure 5. The data 
has been fit with a periodic Yukawa function, 

V(R)=a[ [ - >- + - 'J, (21) 

where TV is the spatial length of the lattice. The parameter d = d(N) pro- 
vides the correct limiting behaviour for r = 0. This form reduces to the con- 
tinuum Yukawa function as N — > oo (d — > 0), and is in excellent agreement 
with numerical calculations of the discrete Fourier transform of fe2 ^ a . We 
found d~0.5atiV = 8, and using this in we obtain (goodness-of-fit 
X 2 /(degree of freedom) = 0.5) 

a = (-12 ±3) x 10~ 2 m Af , 

fx = 0.7 ± 0.4 (lattice units). (22) 

The strength of this potential, described by the constant a, is consistent with 
nuclear potentials which are two orders of magnitude smaller than nucleon rest- 
masses. 

The potential is certainly attractive for small separations. Evidence for 
attraction at small separations was also seen in the staggered-quark calculation 
of Mihaly et. al. [|J . To ascertain whether the potential is significantly non-zero 
for meson spacings R > 0, the data for R = 2 to R = 6 was fit with a constant 
function, 

V(R) = c, (23) 

giving 

c= (-1.3 ±0.8) x 10" 2 m M - (24) 

This constant term is more than one and a half standard deviations below 
zero, suggesting that there is some attractive potential even at medium-range 
separations. We regard the near- vanishing potential at the largest separations 
(R = 4) to be circumstantial evidence that our lattice volume is large enough 
to contain the meson-meson system. 

The Yukawa parameters ( p2| ) correspond to a model of hadron interactions 
mediated by pion-exchange, with the exchanged pion mass given by m, = /i = 
0.7 ± 0.4 in lattice units. Even though we are working in the quenched approx- 
imation and virtual gg-pairs are absent, the exchange term in (JsJ) will include 
pion exchange contributions — one such exchange is shown in Figure 6. 

Pion propagators were constructed using the same gauge field configurations 
at the same value of k, to provide a pion mass for direct comparison with (p2f). 
The pion propagator is equivalent to (||), with two light-quark fields, 

G„(t) =Tr^G 2 (M;0,0), (25) 

n 

and the pion mass was extracted from cosh fits to this propagator. The result 
of TOtt = 0.38 ± 0.18 agrees within errors with the Yukawa pion mass above. 
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4.2 The MM-system 

Simulations of a meson-antimeson system using (^) were performed. The re- 
sults were not consistent with those described above for the MM-system. The 
propagators tended to behave exactly like Wilson Loops, just indicating that 
the system's ground state was a state where the light quarks have annihilated, 
leaving a static QQ pair. Figure 3(b) shows a contribution from such a state. 
A proper simulation of MM systems may require a more realistic, dynamic 
treatment of the heavy quarks. 

5 Conclusions 

The binding potential between two heavy-light mesons in SU(2) QCD was found 
by considering the long-time behaviour of propagators for the two- meson system. 
The heavy-quark approximation simplified analysis considerably, as the heavy 
quark plays a role analogous to the heavy nucleus in the adiabatic calculation 
of the binding potential in the H2 molecule. 

The binding potential derived in this study is attractive and short ranged, 
and periodic due to the boundary conditions of the lattice. The meson-meson 
interaction can be modeled by a pion-exchange Yukawa theory. Comparing 
our results with a Yukawa model, the potential has the correct form, and the 
exchanged pion mass derived from the potential agrees within errors with the 
pion mass calculated at the same values for the simulation parameters. While 
the large errors prompt caution, these results make us confident that the heavy- 
light meson-pair system is tractable in lattice QCD, and is a promising step 
towards a first principles treatment of the nuclear physics problem. 

Unfortunately, in the MM-propagator we only saw evidence of a static QQ- 
pair, due to the annihilation of the light-quarks. Overlap with physical MM- 
states would be enhanced if the heavy-quarks were allowed to propagate spa- 
tially, perhaps in the regime of non-relativistic QCD. More work needs to be 
done to fully appreciate the subtleties of this system. 

The very low light-quark mass used in this study (k = 0.165 gives a rho-pion 
mass ratio of ~ 3.0) resulted in large statistical fluctuations in the quark 
propagators, evidenced by the error bars in Figure 5. Unfortunately, raising the 
quark mass leads to smaller contributions from light-quark exchange, and so 
weakens the resulting potential. These problems are in part due to the low-order 
approximations in the Wilson quark and gluon actions. We have commenced 
simulations using improved actions for the gauge fields ^ ^| and fermion fields 
[ p"ij[ . These actions remove most of the discretisation errors, and restore much 
of the rotational symmetry of the continuum QCD action lost in the hypercubic 
lattice approximation. 
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Q q 

Figure 1: Heavy- light meson propagator constructed from heavy (Q) and light 
(q) quark propagators. 
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Figure 2: (a) Direct and (b) Exchange terms for MM-system. 
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Figure 3: (a) Direct and (b) Exchange terms for MM-system. 
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Figure 4: Fits to the propagator ratio (19), for meson separation R = (circles) 
and R = 1 (triangles). Error bars indicate Jacknife errors. 
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Q q q Q 

Figure 6: A pion-exchange contribution to the MM-propagator — at the time- 
slice indicated by the dotted line, two heavy-light mesons and a pion are present. 
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